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Abstract. We study the finite-size corrections of the dimer model on 00 x N square lattice 
with two different boundary conditions: free and periodic. We find that the finite-size 
corrections depend in a crucial way on the parity of N, and show that, because of certain 
non-local features present in the model, a change of parity of N induces a change of boundary 
condition. Taking a careful account of this, these unusual finite-size behaviours can be fully 
explained in the framework of the c = — 2 logarithmic conformal field theory. 
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1 Introduction 

The dimer model is extremely simple to define. We take a rectangular grid C with M rows 
and N columns, and consider all arrangements of dimers (dominoes) so that all sites of C are 
covered by exactly one dimer. An example of such an arrangement for a 9 x 18 grid is shown in 
Fig. la. In case both M and N are odd, one site has to be removed from the grid; this site is 
conventionally to be one of the four corners. 

One can be more general and introduce monomers (vacancies), namely fixed sites which 
cannot be covered by dimers. A dimer arrangement is possible in presence of m monomers 
provided the numbers of sites which must be covered by dimers and belonging to the even 
sublattice and to the odd sublattice (sum of coordinates even or odd resp.) are equal. In 
particular, MN — m must be even. A dimer configuration on a 9 x 18 grid with four monomers 
is shown in Fig. lb. The case m = (or m = 1 if MN is odd) is referred to as the close-packed 
limit of the dimer model. Corresponding to these arrangements of monomers and dimers, we 
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Figure 1. (a) The left figure shows one among some 4.653 x 10 18 close-packed dimer coverings on a 
9x18 grid, (b) The right figure shows a dimer configuration on the same grid but with four monomers, 
represented by red dots. The insertion of the four monomers has reduced the number of configurations 
by a factor roughly equal to 57.5, down to 8.097 x 10 16 . 



introduce the partition function 

Z{x,y\ Zl ,...,z m )= Yl xUh y nv - (!) 

coverings 

It counts the number of dimer coverings in presence of m monomers located at positions 
zi, . . . , z m , in the bulk or on boundaries, with respective weights x and y assigned to hori- 
zontal and vertical dimers. As the number + n v of dimers is fixed, the partition function 
essentially depends on x, y through the ratio x/y only 

There are obvious generalizations (other types of lattices, coverings by other lattice animals), 
but the basic questions remain the same. Depending on the values of the parameters (here, x/y), 
can the system become critical and does it exhibit phase transitions? What is the probability 
that certain chosen bonds be covered by dimers? How do monomers affect the number of dimer 
coverings? More importantly for our purpose here, what are the finite-size correction terms 
in the partition function, and, in case the system is critical, is its scaling regime conformally 
invariant? 

The dimer model has been originally introduced to model physical adsorption of diatomic 
molecules on crystal surfaces pQ. The first studies of the dimer model were conducted in the 
early sixties, with pioneering works by Kasteleyn [21 [3], Fisher [3], and Fisher and Temperley [5]. 
Soon after that, correlations between dimers and monomers on the square lattice have been 
studied in [HIT], and have been revisited more recently for the triangular lattice in [Hl|9j. The 
effects caused by the insertion of monomers have also been reconsidered recently [101 [TT] I12j . 
Finite-size effects also have a long history, starting in [4j [13], with many subsequent works 

Many critical systems have been shown to have a local scale invariance, so that their scal- 
ing limit can be described by a conformal field theory. Such a theory is primarily (but non- 
exclusively) parametrized by the value of its central charge c, which itself is related to the 
finite-size corrections to the critical free energy. 

The calculation of the central charge based on the finite-size corrections has led to some con- 
fusion in the literature. Indeed it is well-known since [13] that the finite-size corrections depend 
on the parities of M and N. This has prompted several authors to claim that correspondingly 
the central charge was also dependent of the parities of M and iV [lOj . 

Here we reanalyze this question. The finite-size corrections determine univoquely the effective 
central charge c e fj = c — 24/i m ; n [22] . To extract the central charge itself, one needs to determine 
the groundstate energy h m \ n of the Hamiltonian with prescribed boundary conditions. So because 
of /imin the effective central charge depends on the boundary conditions. We show, by changing 
variables from dimer coverings to spanning trees (and more generally, to arrow configurations), 
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that a change of parity of M or N has precisely the effect of changing the boundary conditions. 
The known values of /i m i n and the values of c c g as computed from the finite-size corrections 
enable us to obtain a consistent value for the central charge, here equal to c = —2. Surprisingly, 
this is not the only consistent value of c that can be used to describe the dimer model, since it 
has been shown in that a dimer configuration can be encoded in a height function [23], which 
converges in the scaling limit to a Gaussian field, with c = 1 [24]. We will discuss this peculiar 
situation in the Conclusion. 

In what follows, we consider the dimer model on a square lattice in the close-packed limit, 
and for two different boundary conditions, free and periodic. In Section 2, we collect the known 
results for the finite-size corrections, and the way these relate to the central charge. In Sections 3 
and 4, we compute the central charge from the finite-size analysis respectively on a strip with 
free boundary conditions, and on a cylinder (periodic boundary conditions). The last section 
summarizes our point of view on the description in terms of c = —2 versus the one based on 
c = 1, relating this issue to the more general monomer-dimer problem. The present results have 
been reported in [21], but we take here the opportunity to give more details. 



2 Finite-size analysis 

The partition function ([1]), in the close-packed case, has been first computed by Kasteleyn as 
the Pfaffian of a certain matrix [2]. Since then it has reproduced by a variety of methods [25] . 
including that based on spanning trees |26[ [27] which we will use later on. 

The partition function, up to an irrelevant factor, depends on x/y only. In the infinite volume 
limit, the dimer model on the square lattice becomes critical but does not show a phase transition 
(unlike on other lattices). The critical properties do not depend on the values of x/y (provided 
it is neither zero nor infinite), so we set x = y = 1 in the following. The partition function, 

Z(M,N) = # dimer coverings, 

simply counts the number of ways the grid C (minus a corner if MN is odd) can be fully covered 
by dimers. The topology of C is fixed by the boundary conditions: it forms a rectangle if free 
boundary conditions are imposed in two directions (like in Fig. la), or a cylinder if periodic 
boundary condition is chosen in the horizontal direction. 

We are especially interested in the free energy Fn per unit of height for anc»xiV lattice, 

F N = - lim -J- log Z(M,N). 

M— >oo 1V1 

The free energy depends on the boundary conditions, free or periodic, and on the parity of N. 
This leads to four quantities F^ e e ven , F^ e dd , F^ vcn and F^ odd . 

The partition functions with these boundary conditions can be expressed in terms of Z a g for 
a,/3 = 0,| [20], with 

N-l AI-l 

zz J} {i,M,N)= n n 4 

n=0 m=0 

The full expansion of Z a g for large M, N has been obtained in [19]. Using these results, 
one can easily obtain the dominant finite-size corrections terms for the above four free energy 
densities, as power series in N. 

For an infinitely long strip of width N with free boundary conditions, one obtains 

^ = ~W-f + ilo g (l + ^ + ^i + ..., (2) 



. 2 vr(n + a) . 2 vr(m + 0) 
sin — h sin 



N 



M 
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F^ en = --N-^ + hog(l + V2)-^ + ... , (3) 

7T TT Z Z4 iV 

where G = 0.915965 is the Catalan constant. 

The analogous results for the periodic case, i.e. an infinite cylinder of perimeter N, read [13j 

^npr G „ T 7T 1 „npr G TT 1 , . , 



The dominant term in these expressions yields the exponential growth of the number of dimer 
coverings of a M x N grid, namely Z(M,N) ~ ^P/^MN = [i,jgiQ2) MN , for both types of 
boundary conditions. 

The general form of the free energy per unit length of an infinitely long strip of finite width N 
at criticality is 

F N = /buikiV + 2/ surf + - + ..., (5) 

where /bulk is the free energy per bulk site, and /bulk + /surf is the free energy per boundary 
site (assuming identical boundary conditions on the two sides of the strip), and A is a constant. 
Unlike the free energy densities /bulk an d / SU rf ; the constant A is universal. The value of A is 
related to the central charge c of the underlying conformal theory, and depends on the boundary 
conditions in the transversal direction. The result is that A is proportional to the effective central 
charge c cff = c - 24/i min [28l l29l 130] . 

A = c efr = 7r (/i min - -^j on a strip, (6) 

7T / C \ 

A = - - c eff = 47r [h min - —J on a cylinder. (7) 

The number /i mm is the smallest conformal weight in the spectrum of the Hamiltonian with the 
given boundary conditions (for the cylinder, we assumed that this operator is scalar, h m \ n = 

hmin) • 

The values of A given above are most easily understood from the transfer matrix formalism 
for calculating the lattice partition function, even though the result seems generally more valid 
(i.e. when there is no transfer matrix). Let us first consider a cylinder of perimeter N and 
height M, with certain boundary conditions \b) and \t) on the bottom and top edge ail. The 
transfer matrix T cy \ is labelled by the degrees of freedom living on a horizontal circle, in terms of 
which the partition function is equal to Z(M,N) = (b\T c ^\t). In the thermodynamic limit, the 
vector space spanned by the degrees of freedom living on a circle (space coordinate) goes over to 
an infinite Hilbert space 7i. cy \, and the transfer matrix, which can be viewed as a unit translation 
operator in the vertical direction (time coordinate), can be written as T cy \ = e~ Hc y l , in terms of 
a Hamiltonian H cy \. In the large M limit, the partition function Z(M,N) = (6|e _M ^ cyl |t) will 
be dominated by ground-state Eq of H cy \, so that F N = - lim fa log Z(M, N) = E . 

The Hamiltonian is the charge associated to the time-time component of the stress-energy 
tensor H cy \ = ^ duToo- If the system is critical and conformally invariant, the Fourier modes 
of the stress-energy tensor are the left and right Virasoro modes L n , L n , and the Hamiltonian 
is simply given by the zero-th moded Virasoro generators as H cy \ = ^f(Lo + Lq — ^). This last 
formula assumes a normalization where the ground-state energy vanishes in the thermodynamic 
limit N — ► oo, and therefore ignores the bulk term in ([5]) (in the periodic geometry we consider 



1 In the limit M — > oo, we expect that the specific boundary conditions we choose do not matter if the theory 
is local, since the boundaries are sent off to infinity. We could as well choose periodic boundary condition in the 
vertical direction, closing the cylinder into a torus. 
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here, there is no surface and therefore no surface term in ([5])). In addition the Hilbert space 
decomposes into representations of the left and right Virasoro algebras as TC cy i = (B h ^ N h ^ TZh <8) 
IZfr, in which the N h ^ are integer multiplicities. The representations TZh<S>7Zh are highest weight 
representations, meaning that they are built from a highest weight state \h) ® \h) by applying 
all Virasoro modes. All states of the representation are eigenvectors of Lq + Lq, but the state 
with the smallest eigenvalue is the highest weight state, (Lq + Lq — h — h)\h) ®\h) = 0. Putting 
all together, one obtains that the ground-state of H cy \ in TC cy i is equal to 

_ 27T / 7 C \ 

Lq — -jj: ^/imin T" "-min ^0/ ' 

where h m [ n and h m { n label the representation in T~C cy i with the smallest (Lq + Lq) eigenvalue. 
Assuming the equality h m \ n = h m \ n yields the result quoted above. 

For the strip, a first change is that one has to specify the boundary conditions a, b on the 
edges on the strip. The lattice transfer matrix depends on a, b, as do the Hilbert space 7~i a ,b an d 
the Hamitonian H a i, in the thermodynamic limit, but the formula Fn = E^ ,b still holds. 

The second change is that there is only one copy of the Virasoro algebra on the strip. It implies 
that the Hamiltonian is now equal to H a i, = j^(Lq — ^), and the Hilbert space decomposition 

reads H a> b = ©h-^V?' IZh, leading directly to the value of A given in ([6]). Again the normalization 
means that the quantum Hamiltonian ignores the bulk and surface terms in ([5]). 

The finite-size corrections computed above for _F free and F per have the correct form ([5]), but 
at first sight look paradoxical. Comparing them with Q and (J7J), we see that the effective 
central charge depends on the parity of N. The boundary conditions do not seem to change 
with the parity of N. This would imply that h m [ n does not change either, that the central charge 
itself has to change. This would be most peculiar since the central charge also controls the bulk 
conformal theory. 

The only way out is to accept that the boundary conditions change with the parity of N, 
although this is not apparent in the dimer variables. We will show in the following that the 
effective central charge changes with the parity of N, c c g = —2 for N odd and c e g- = 1 for N 
even, not because the central charge changes, but because the value of /i m j n changes, due to the 
fact that a change in the parity of N effectively changes the boundary condition. On the strip, 
this effect has been already noted in |31| . 

To understand this peculiarity of the dimer model, we consider a change of variables, namely 
we replace dimer configurations by arrow configurations on a sublattice. On the strip, the 
arrow configurations define spanning trees, so that the dimer model is mapped to the spanning 
tree model [271 EI] or, equivalently, the Abelian sandpile model [32]. The case of the cylinder 
is slightly more complicated because the arrow configurations do not always define spanning 
trees. The analysis of this case is however similar. 

3 Dimers on a strip 

Let us consider first the dimer model on the rectangular lattice C of size M x N with free 
boundary conditions. Since we are interested in the limit M — > oo, the parity of M will not 
matter here. For simplicity, we take M odd, and discuss successively the cases N odd and N 
even. 

When M and N are both odd, there is a well-known bijection [26] between close-packed dimer 
coverings of C with one corner removed and spanning trees on the odd-odd sublattice Q C C 
(thus Q contains the sites whose coordinates, counted from the lower left corner, are both odd). 

A dimer containing a site of Q, in red in Fig. 2, can be represented as an arrow directed along 
the dimer from this site to a nearest neighbour site in Q. It is easy to prove that the resulting 
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Figure 2. The left figure represents a dimer covering of a 9 x 17 lattice. The red dimers are those which 
touch the odd-odd sublattice Q. The right figure shows the corresponding spanning tree on the 5x9 
lattice G, with a unique root pictured by an open dot. 



set of arrows generates a uniquely defined spanning tree; all the arrows point to a unique root, 
located at the corner which had been removed from C (see Fig. 2). Since the dimers which do 
not contain a site of Q are completely fixed by the others, one has a one-to-one correspondence 
between dimer coverings on C minus a corner and spanning trees on Q, rooted at the removed 
corner. The Kirchhoff theorem then expresses the number of dimer configurations as Z = det A, 
where A is the Laplacian matrix on Q with appropriate boundary conditions, see Section IV for 
a proof of this result. Viewing Q as a graph, one sets Ajj = — 1 if sites i and j are connected 
in Q, and An is equal to the number of sites in Q that i is connected to, plus 1 if i is connected 
to the root. As shown in [32], spanning trees on Q, rooted at a corner, are in bijection with the 
configurations of the Abelian sandpile model (ASM) on Q, with closed boundary conditions on 
the four boundaries, the only sink (dissipative) site being the root of the trees. 

When M — > oo, the lattice Q becomes an infinitely long strip of width N. The root is sent 
to infinity so that no boundary site (at finite distance) is connected to the root. In the ASM 
language, this means that the boundaries are not dissipative and hence there is no out-flow of 
sand at the boundaries; such boundaries conditions have been called closed boundary conditions. 

There are now very strong arguments to believe that the (scaling limit of the) ASM on 
a square lattice is described by a logarithmic conformal field theory with a central charge c = — 2 
[321 l33l IMj [35l 136] . In particular, the spectrum of the ASM Hamiltonian on a slice of the strip 
with closed boundary conditions at the two ends has been computed in [33]. There are two 
ground-states, the identity operator and its logarithmic partner, both of conformal weight 0, so 
that /i m in = 0. The effective central charge in this sector is therefore c e g = —2, and the general 
formula ([6]) reproduces the finite-size corrections (|2|). 

When M is odd and N is even, dimer coverings exist without the need to remove a corner. 
In this case, the above construction leads to a set of spanning trees on the odd-odd sublattice Q, 
where certain arrows may point out of the lattice from the right vertical side, see Fig. 3. Viewing 
this vertical boundary of Q as roots for the spanning trees, we see that dimer coverings on C 
map onto spanning trees on Q which can grow from any site of the right side. The sites on this 
boundary, being connected to roots, are dissipative in the ASM language, and form an open 
boundary. Thus the spanning trees map onto the ASM configurations with one vertical open, 
dissipative boundary, and the three other closed. 

In the limit M — > oo, the lattice becomes an infinite strip with open and closed boundary 
conditions on the two sides. In this case, the ground-state of the Hamiltonian with such boundary 
conditions is a primary field of conformal weight /i m i n = — 1/8 [33]. With c = — 2, this yields 
c e ff = 1 and again the formula ([6]) gives the correct result ([3]). 

Thus the leading finite-size corrections for an infinitely long strip of width N agree with the 
prediction of a c = —2 conformal field theory, provided one realizes that changing the parity 
of N genuinely changes the boundary conditions, an effect due to the strong non-locality of the 
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Figure 3. The same dimer covering of a 9 x 18 lattice as in Fig. 1 is shown on the left, where the dimers 
which touch the odd-odd sublattice Q are coloured in red. The right figure shows the corresponding 
spanning tree on the 5x9 lattice Q ; the sites on the right boundary are all connected to roots (the open 
dots), although in this particular instance, two of these connections have not been used by the tree. 



O 



Q 



Figure 4. (a) The left figure shows the spanning tree associated to the dimer covering of Fig. 3, drawn 
on the even-even sublattice rather than on the odd-odd sublattice (the arrows are no longer shown). 
All boundaries are open (their sites are connected to roots) except the right boundary, which is closed, 
(b) The same spanning tree is reproduced on the right figure, where in addition the spanning tree drawn 
on the odd-odd sublattice, in blue, is superimposed. The duality of the two spanning trees is then 
manifest. 



dimer model [21]. The change of boundary conditions is not apparent in the dimer model itself, 
but is manifest when one maps it onto the spanning tree model or the sandpile model. 

Interestingly the bijection between the dimer coverings and the spanning trees holds if we 
use the even-even, even-odd or odd-even sublattice. The boundary conditions however change, 
but one can easily see what the changes are. If, instead of choosing the odd-odd sublattice, one 
selects the sites whose horizontal (resp. vertical) coordinates are even, the left and right (resp. 
bottom and top) boundaries change from closed to open, and from open to closed. So if N 
is odd, the vertical sides become open rather than closed. The spectrum of the corresponding 
Hamiltonian changes, with a non-degenerate ground-state being the identity operator [33], so 
that the value h m \ n = remains. If N is even, the left and right boundaries, previously closed 
and open respectively, become open and closed respectively, so the Hamiltonian remains the 
same, h m [ n = —1/8. Fig. 4 shows the spanning tree associated to the dimer covering in Fig. 3 if 
we choose the even-even sublattice. 

In fact the odd-odd and even-even sublattices are dual, as are the corresponding spanning 
trees, see Fig. 4b. The same is true of the odd-even and even-odd sublattices. Under this duality, 
the boundary conditions open and closed are exchanged. 

4 Dimers on a cylinder 

We consider here an M X N rectangular lattice C with periodic boundary condition in the 
horizontal direction, so that £ is a cylinder of perimeter N and height M. As before, we will 
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eventually take M to infinity, which makes its parity irrelevant. For convenience we choose M 
even. We discuss the cases N odd and iV even separately. 

If N is odd, we select again the odd-odd sublattice Q. It is easy to see that two columns 
of Q will be neighbours in Q and in C (connected by horizontal bonds). Therefore a dimer may 
cover zero, one or two sites of Q. As before the dimers covering no site of Q are completely fixed 
by the others and play no role. For the others, we do the same construction as in the previous 
section. If a dimer touches only one site of Q, we draw an arrow directed along the dimer from 
that site to the nearest neighbouring site of Q. However, for a dimer laid on two sites of Q, the 
two arrows would point from either site to the other, ruining the spanning tree picture. It can 
nevertheless be restored in the following way. 

Instead of seeing the two arrows as pointing from one site to its neighbour, we say that they 

point towards roots inserted between the neighbour sites, thus replacing the arrows •-» «— • 

by »-K3 o+-%- This in effect amounts to opening the cylinder by removing the horizontal bonds 
of C which connect sites of Q, unwrapping it into a strip, and to adding columns of roots on 
the left and on the right side of the strip. The new arrow configurations define spanning trees, 
rooted anywhere on the left and right boundaries. So dimer coverings on the original cylinder 
are mapped to spanning trees on a strip, with open upper and closed lower horizontal boundaries 
(by the arguments of the previous section, since we chose M even), and open vertical boundaries. 

However these two vertical boundaries are not independent. If a tree grows from a site i on 
the left boundary, then a tree must grow from the site j on the right boundary, where j is the 
site which were neighbour of i in Q before the opening of the cylinder. Conversely if no tree 
grows from i, then no tree can grow from j. Thus the growth patterns on the left and on the 
right boundaries must exactly match. 

For finite N, this correlation between the two boundaries will be felt throughout the lattice. 
As N (and M) increase, the generic properties of the spanning trees will not change if the trees 
are allowed to grow from slightly different boundary sites. Thus if the density of roots on the 
left and on the right boundaries is finite, it will make no difference whether the tree growth 
patterns are identical or random on each boundary. 

Following this argument, when M goes to infinity, the lattice becomes an infinite strip 
with open boundary condition on either side. As mentioned above, the ground-state of the 
Hamiltonian is the identity, of weight h m { n = 0, leading to an effective central charge c e g = — 2. 
The general formula (jfij) for the strip gives the correct result (j3J). 

This is a very unusual situation. Although the dimer model is originally defined on a cylinder, 
it shows the finite-size corrections expected on a strip, and must really be viewed as a model on 
a strip. 

For N even, the problem of having a dimer occupying two neighbouring sites of Q does not 
arise. Hence the arrows point from one site of Q to a neighbouring site of Q and they never 
overlap. However the arrow configurations one obtains do not define spanning trees because 
there can be sequences of arrows looping around the cylinder. In general the arrows form 
a combination of oriented loops wrapped around the cylinder with trees growing from the loops. 
The one-to-one correspondence between the oriented loops combined with tree branches on one 
side, and dimer configurations on the other side can be established as above. The enumeration 
of the loop-tree configurations requires a slight generalization of Kirchhoff 's theorem. In order 
to appreciate it, we make a little digression, to recall the combinatorial content of the usual 
Kirchhoff theorem (see for instance [21])) as use d in Section 3. 

Define the Laplacian of an unoriented graph Q as the symmetric matrix A with entries 
Ajj = — 1 if there is a bond connecting site i and site j ^ i, and Z{ = Aj i equal to the number 
of sites to which % is connected, including the root(s). Thus ^ • Ajj is the number of connections 
from i to the root(s). Let us now place, at each site i of Q, an arrow pointing to any one of the 
Zi neighbours of i (possibly to the root if i has a connection to it). Such a set of arrows defines 
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an arrow configuration on Q, and we want to show that the determinant of A precisely counts 
those arrow configurations which contain no loop. 

The determinant is a sum over the permutations in SV with iV = the number of sites 

of g, 

det A = ^ £a A 1,<t(1) A 2,a(2) ' ' ' &N,<t(N) ■ 
o-gSjv 

Writing a permutation as a product of cycles of lengths £ m , and using e a = (— l) Ar +# c y cles ) it is 
not difficult to see that each term in the sum carries a sign equal to (— l)#P r °P cr c y clcs where by 
proper cycle we mean a cycle of length strictly larger than 1. Up to this sign, a permutation a 
with exactly k proper cycles brings a contribution equal to F] i z% where the product is over those 
sites left invariant by the permutation. Combinatorially this number counts the arrangements 
of arrows containing k loops whose locations in the graph are completely specified by the cycle 
structure of a; on the other hand, the arrows coming out of the sites fixed by a are free to point 
in any of the allowed directions, and may therefore form themselves other loops. 

The sum over the permutations can now be reorganized as a sum over the number k of proper 
cycles, 

[iV/2] 

detA=^(-l) fc Yl \A lia{1 yA NMN) \. (8) 

k=0 a has k proper cycles 

The term for k = 0, equal to Ylieg Zi coun ts the total number of unconstrained arrow 
configurations. The term k = 1, up to the minus sign, counts the arrow configurations which 
have at least one loop of a fixed type, and sums over the possible loops. It therefore overcounts 
the number of arrow configurations with at least one loop by the number of arrow configurations 
which have at least two loops. This is taken care of by the k = 2 term, which however overcounts 
it by a amount equal to the number of configurations with at least three loops, and so on. By 
the inclusion-exclusion principle, the alternating sum ([8]) exactly counts the configurations of 
arrows with no loop at all. 

Let us now come back to the problem of counting the loop-tree configurations of arrows, 
in which the only loops we allow must wrap around the cylinder (i.e. be non-contractible) . 
Because they cannot intersect themselves, the loops can wrap only once around the cylinder. 
For this, we modify the previous Laplacian A by changing from —1 to +1 the entries Ajj for all 
pairs i,j such that the bond i — j crosses a line going from one boundary of the cylinder to the 
other boundary, which we will call a defect line. The resulting still symmetric matrix A^ can 
be viewed as the antiperiodic Laplacian (the lower and upper boundaries are closed and open, 
or vice- versa, since we have taken M even) . We show that Z = det A a precisely counts the 
loop-tree configurations we need. 

The key observation is that a contractible loop on the cylinder crosses the defect line an 
even number of times, whereas a non-contractible loop crosses it an odd number of times. As 
a consequence, and with respect to the previous situation, there is an extra minus sign for each 
proper cycle giving rise to a non-contractible loop of arrows. One thus obtain 

[N/2] 

detA A = ]T(-l) fc Yl (-l) #non - contr -|A liCT(1) ..-A w) |. 

k=0 a has k proper cycles 

If we write k = p+q, where p is the number of non-contractible loops (NCL) and q is the number 
of contractible loops (CL), the summation over k is replaced by two summations over p and q, 

detA A = Y/Z(-^ q Yl \^i,a(i) ■ ■ ■ A NMN) \. 

p>0 q>0 a has p NCL.q CL 
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One sees that the arrow configurations with a fixed number p of non-contractible loops, i.e. the 
terms with q = 0, are all counted positively. Moreover, for each p, the inner summation over q 
is an alternating sum which implements the inclusion-exclusion principle and removes all con- 
tractible loops. Therefore det Aa exactly counts the arrow configurations with no contractible 
loops, as claimed. 

In the continuum limit, it becomes the partition function of a free theory of antiperiodic 
Grassmannian fields which, in turn, gives c = —2 and h m [ n = — 1/8 [33J . This is again consistent, 
since the finite size correction @ together with the general formula ([7]) for the cylinder yield 
Ceflf = 1. 

5 Conclusion 

By analyzing the finite-size effects in terms of the effective central charge c e ff = c — 2Ah ra \ n , 
we have shown that the non-local boundary effects in the close-packed dimer model can be 
consistently accounted for by a single conformal theory having central charge c = —2. We 
have provided a consistent framework for understanding the dependence of the finite-size effects 
upon the boundary conditions. However this should not be taken as a proof that c must be 
equal to —2. Indeed since the effective central charge merely determines some combination of c 
and h m i n , one cannot obtain the values of both without some assumption about one of them. 
This assumption can be a posteriori justified if the conformal description obtained from it is 
fully consistent. 

It turns out in this case that another consistent conformal description exists, with c = 1 
[23\ 124] , although a detailed analysis of boundary conditions and parity dependence effects has 
not been carried out in this context. Our explanation for this curious fact is that the c = 1 
theory not only describes the close-packed dimer model, but the general monomer-dimer model 
(it has been suggested in [37] that trimers would require c = 2). 

On one hand, this has been illustrated for instance in [8] where the monomer 2-point cor- 
relation function [6] has been interpreted in terms of two uncoupled Majorana fermions. The 
same interpretation can be made for the general n-point function for monomers on a bounda- 
ry [HE]- On the other hand, the spanning tree description leading to the value c = —2 cannot 
describe dimers with monomers in generic positions, for the basic reason that they are defined 
on a sublattice and therefore cannot keep track of the positions of all the monomers. 

Thus the conformal theory with c = —2 must be viewed as a subtheory of that with c = 1, 
as it is able to describe the degrees of freedom corresponding to dimer coverings but not general 
monomer insertions. For those degrees of freedom, the two descriptions should be equivalent. 
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